i 


^  i  a 


AD-A230  648  - 

REPORT  DOCUMENTATION  PAGE 


vU!  X 


Form  Approved 
OMB  No.  0704-0188 


Public  reporting  burden  for  this  collection  of  information  Is  estimated  to  average  1  hour  per  response.  Including  the  time  for  reviewing  Instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  Information  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  Information  Including 
suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services.  Directorate  for  Information  Operations  and  Repons.  1215  Jefferson  Davis  Highway.  Suite  1204.  Arlington.  VA  22202-4302. 
ana  to  the  Office  of  Management  and  Budget.  Paperwork  Reduction  Project  (0704-0188)  Washington.  DC  20503 


3  REPORT  TYPE  AND  DATES  COVERED 

Professional  paper 

5  FUNDING  NUMBERS 

PE:  ZW18 
WU:  DN309  170 
PE:  0G01152N 


PERFORMING  ORGANIZATION 
REPORT  NUMBER 


1  AGENCY  USE  ONLY  (leave  blank) 

2  REPORT  DATE 

December  1990 

4  TITLE  AND  SUBTITLE 

LIQUID  METAL  FUEL  COMBUSTION  MECHANICS 

6  AUTHOR(S) 

T.  Duvvuri 

7  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Naval  Ocean  Systems  Center 

San  Diego,  CA  92152-5000 

9  SPONSORING/MONITORING  AGENCY  NAME|S)  AND  ADDRESS(ES) 

Office  of  Chief  of  Naval  Research 

Independent  Research  Programs  (IR) 

OCNR-lOP 

Arlington,  VA  22217-5000 

10  SPONSOR ING/MONITORING 
AGENCY  REPORT  NUMBER 


12a  DISTRIBUTION/AVAILABILITY  STATEMENT 


12b  DISTRIBUTION  CODE 


Approved  for  public  release;  distribution  is  unlimited. 


13  ABSTRACT  (Maximum  200  words) 

\  f 

The  modeling  of  the  droplet  formation  at  the  gas/liquid  boundary  interface  of  a  gaseous  jet  injected  into  a  liquid  metal  bath 
and  the  turbulent  mixing  of  the  resultant  two-phase  (gas/liquid  droplets)  mixture  is  presented  as  a  preliminary  to  the  analysis  of 
the  liquid  metal  fuel  combustion  problem.  The  model  is  used  to  predict  velocity  and  liquid  droplet  fraction  distributions  across  the 
mixing  zone  and  along  the  centerline  of  the  jet.  These  results  show  that  the  distribution  of  the  liquid  fuel  droplets  peaks  away 
from  tlie  region  of  maximum  oxidizer  concentrations.  ^  -j* 


DTIC 

i  ELECTEI 
JAN  10 1991)1 

1  B  I 


* 

D 


Published:  A1AA/SAE/ASME/ASEE ,  26th  Joint  Propulsion  Conference,  July  1990 


17  SECURITY  CLASSIFICATION 
OF  REPORT 

UNCLASSIFIED 


18  SECURITY  CLASSIFICATION 
OF  THIS  PAGE 

UNCLASSIFIED 


19  SECURITY  CLASSIFICATION 
OF  ABSTRACT 

UNCLASSIFIED 


15  NUMBER  OF  PAGES 


20  LIMITATION  OF  ABSTRACT 


SAME  AS  REPORT 


NSN  7540-01-280  5500 


Standard  form  298 


At  A  A 


AIAA  90-2308 

LIQUID  METAL  FUEL  COMBUSTION  MECHANICS 
TIRUMALESA  DUVVURI 
Naval  Ocean  Systems  -center 
San  Diego,  CA 


|  Aooesslon  For  .  ) 

NTIS  GRAAI 

DTIC  TAB 

□ 

Unannounced 

□ 

Ju6tlf leal ion _ 

By - 

Distribution/ _ 

Availability  Codes 
jA^all  and/or 
Dlst  I  Spoolal 


AIAA/SAE/ASME/ASEE 
26th  Joint  Propulsion  Conference 

July  16-18,  1990 /Orlando,  FL 


*  V 


For  permission  to  copy  or  republish,  contact  the  American  Institute  of  Aeronautics  and  Astronautics 
370  L’Enfant  Promenade,  S.W.,  Washington,  D.C.  20024 


oV 


LIQUID  METAL  FUEL  COMBUSTION  MECHANICS 


Tirumalesa  Duwuri* 
Naval  Ocean  Systems  Center 
San  Diego,  CA  92152-5000 


Abstract  jo  Density 

' 6*  Viscous  stress  tensor;  surface  tension 

The  modelling  of  the  droplet  formation  at  the  gas/liquid  Turbulence  scale  for  droplet  dispersion 

boundary  interface  of  a  gaseous  jet  injected  into  a  liquid  metal  ^  Differential  Operator,  grad 
bath  and  the  turbulent  mixing  of  the  resultant  two-phase 
(gas/liquid  droplets)  mixture  is  presented  as  a  preliminary  to  Subscripts 
the  analysis  of  the  liquid  metal  fuel  combustion  problem.  The 

model  is  used  to  predict  velocity  and  liquid  droplet  fraction  d  Liquid  Droplets 

distributions  across  the  mixing  zone  and  along  the  centerline  of  g  Gas 

the  jet.  These  results  show  that  the  distribution  of  the  liquid  k  Species  no.,  droplet  no. 
fuel  droplets  peaks  away  from  the  region  of  maximum  oxidizer  1  Liquid 

concentrations.  r  Relative 


nomenclature 

a  Speed  of  sound 

Ba  Empirical  constant  for  atomization  drop  size 
Ca  Empirical  constant  for  atomization  rate 
Cy  Specific  heat  at  constant  volume 
d  Liquid  Droplet  diameter 

Droplet  drag  per  unit  relative  velocity 
F  Fractional  volume  of  gas 
G  Body  force 
h  Specific  Enthalpy 
H  Liquid  Droplet  specific  enthalpy 
I  Specific  internal  energy 
K  Thermal  conductivity 

Kp  Turbulence  model  empirical  constant,  =  0(0.1) 
m  Droplet  mass 

Mg  Mass  stripped  away  at  the  gas/liquid  interface 
r7T'  Molecular  weight 
Mu  Mussel t  numbert 
p  Pressure 
Pr  Prandtl  nutter 

r  Droplet  radius;  or  radial  coordinate 
R  Universal  gas  constant 
Re  Reynolds  number 
S'  Droplet  momentum  coupling  term 
Sc  Schmidt  rtunber 

SH  Droplet  angular  momentum  coupling  term 
t  Time 
T  Temperature 

u,v,w  Gas  velocity  components 

U, V,W  Liquid  droplet  velocity  components 

Vm  RMS  vetocity  of  gas  molecules 
We  Weber  number 

x  Radial  or  transverse  coordinate 
y  Axial  coordinate 

£  Turbulent  kinetic  energy  level  for  droplet  dispersion 
A  Cell  scale  length  for  turbulence  model 
u  Viscosity 


INTRODUCTION 

The  impetus  for  Liquid  Metal  Fuel  Combustion  research  arose  from 
an  interest  in  finding  closed  loop  high  energy  density  compact 
propulsors  for  underwater  applications.  A  Stored  Chemical  Energy 
Propulsion  System  (SCEPS)  produced  by  combining  a  closed-volume, 
liquid  metal  combustor  with  the  boiler  of  a  classical  Rankine  cyc|e^ 
shown  schematically  in  Fig.  1,  satisfies  these  requirements. 
Fig. 2  shows  the  schematic  of  the  boiler-reactor,  and  Figs.  3  and  4 
show  two  types  of  the  metal  combustor  itself. 

In  either  case  a  gaseous  oxidizer  jet  is  injected  (fig. 5)  into 
a  liquid  metal  fuel  bath  in  the  combustor.  As  the  oxidizer  and  the 
fuel  are  separately  introduced,  the  combustion  process  in  this 
case  is  mixing  limited.  The  fundamental  mechanisms  involved  in  the 
combustion  process  for  such  a  case  (gaseous  oxidizer  jet  into 
a  liquid  fuel)  are  not  well  understood. 

Presently  available  studies  analyze  tht  oroblem  under  the 
Locally  Homogoneous  Flow  (LHF)  approximation.  y  These  analyses  do 
not  consider  the  important  processes  that  occur  at  the 
gas/liquid  interface  as  well  as  in  the  mixing  region.  Namely,  (A) 
the  liquid  droplet  formation,  their  agglomeration  and  breakup 
as  well  as,  (B)  the  phase  changes,  (vaporization  and 
condensation),  and  (C)  their  effect  on  the  combustion  process.  (D) 
Further,  the  effects  of  velocity  and  thermal  lags,  (between  the 
gaseous  phase  and  liquid  droplets),  on  the  turbulent  mixing  and 
hence  on  the  efficiency  of  the  (mixing  limited)  ombustion  are 
not  known. 

It  is  well  knowot  from  fuel  spray  combustion  studies  for  gas 
turbine  combustors,  that  fuel  droplet  sizes  can  have  direct 

impact  on  almost  all  key  aspects  of  combustion.  The  initial  droplet 
sizes  formed  at  the  injector  exit  and  their  trajectories,  coupled 
with  the  combustor  primary  zone  turbulent  mixing  significantly 
influence  the  evaporation  and  combustion  characteristics  of  the  fuel 
droplets.  The  droplet  trajectories  themselves  are  influenced  by  the 
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velocity  and  thermal  lags  between  the  phases.  Also  the  droplet 
evaporation,  ignition  and  combustion  is  also  dependant  on  the 
droplet,  size  distribution  (amount  of  surface  available  for 
‘efficient  heat  transfer)  and  its  thermal  lag  (temperature 
differential  between  the  phases)  from  surrounding  gas. 

Previous  studies,  by  the  author,  of  turbulent  mixing, 
ignition  and  combustion  of  two  phase  (gas-solid-particle/ 
pure-gas)  flows  have  shown  these  velocity/thermal  lags  to  play 
an  important  part  in  the  slow  down  of  the  spread  of  the  two 
phase  mixture  (in  particular  the  solid  fuel  particles)  into  the 
mixing  zone1-3  3  and  thus  keeping  the  center  two-phase  jet 
intact  to  a  longer  axial  length.  (This  inefficient  mixing  of 
the  oxidizer  and  the  fuel  leads  to  combustion  inefficiencies  due 
to  non-  optimun  oxidizer/fuel  mixture  ratios.)  Similar  observa¬ 
tions,  about  longer  intact  center  jet,  have  been  made  in  recent 
experimental  studies  at  NOSC  by  Parnell  et  al  .  These  studies 
show  "a  continuous  gaseous  jet  extending  about  60X  of  the 
length  of  the  combustion  chamber",  whereas  calculations  by  LHF 
aproximation  analyses  by  Lin  etal  of  PennstateV  ,  show  a  much 
shorter  jet  core. 

The  aim  of  this  study  is  to  provide  an  insight  into  these 
fundamental  mechanisms.  Understanding  of  these  fundamental 
mechanics  of  the  mixing,  igni t ion  and  combust  ion  process  would 
aid  in  designing  very  efficient  propulsion  devices.  This  would 
also  further  advance  our  understanding  of  the  fundamental  mecha¬ 
nics  of  the  multi -phase  multi -specie  turbulent  reacting/  combus¬ 
ting  flows. 

This  is  a  continuing  research  program  and  this  paper  presents 
the  progress  to  date.  As  a  first  step  in  analyzing  this  complex 
problem,  the  present  work  considers  the  turbulent  mixing  of  the 
two-phase  (gas/liquid  particle)  non- combusting  mixture.  The  work 
accomplished  includes  a  critical  review  of  previous  work, 
modelling  of  the  number  density  distribution  of  the  liquid  metal 
fuel  droplets  generated  at  the  gas/liquid  boundary  as  well  as  the 
velocity  lag  effect  between  the  gaseous  and  liquid  droplet  phase. 
This  model  was  incorporated  into  fluid  flow  equations  for  the 
turbulent  mixing  of  a  two-phase  (gas/liquid  droplet)  mixture.  The 
results  of  some  preliminary  calculations  are  presented  and 
discussed. 
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The  gas  mixture  obeys  the  summation  rules  of  : 
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THE  GOVERNING  EQUATIONS 


There  are  three  fluid  phases  that  need  to  be  considered  in 
our  problem,  namely,  the  gas,  the  liquid  and  the  droplets.  The 
equations  for  the  gas  and  liquid  are  the  usual  Eulerian  continuum 
fluid  dynamic  equations  augmented  by  the  mass,  momentum  and 
energy  exchange  terms  due  to  the  droplet  phase.  The  droplets 
are  treated  as  discrete  particles  travelling  through  the 
continuum  ajjd  their  dynamic  behaviour  is  described  by  Lagrangian 
equations. 


and  the  mass  conservation  equation  for  species  k  is  given  by 
^Pfc  .  U*  V7  v  .  1  P) 


t^.V(rPy 


(8) 


For  our  case  where  a  gaseous  jet  is  injected  into  a  liquid 
bath,  the  liquid  enters  into  the  picture  mainly  at  the  gas/liquid 
interface  where  a  radial  velocity  is  induced  to  the  fluid  due  to 
the  entrainment  effect  of  the  gaseous  jet.  The  liquid  again 
enters  into  the  picture  once  the  gas  jet  has  dissipated  its 
momentum  and  subs  tried  by  the  liquid  or  condenses  to  liquid. 

The  gaseous  phase  will  be  considered  to  be  consisting  of  a 
multi-species  compressible  mixture  whi  le  the  l  iquid  phase  wi  1 1  be 
assumed  to  be  incompressible  with  a  constant  density.  The  mass 
conservation  equation  for  the  two-phase  fluid  is: 


Summing  over  all  species,  the  mass  conservation  equation  for  the 
gaseous  phase  alone  is, 


5 1  r 


V(rPg) 


I  F 


-h(9) 

F 


5 p/dt+-jV.  (rpu)  -ps 


where is  the  mass  of  liquid  that  is  stripped  away  as  droplets 
or  as  a  result  of  phase  change  (evaporation  or  condensation), 
and 


Combining  equations  1,2,9,  the  equation  for  F  is  obtained  as: 


4^  +  -V.  ( ruF )  --V.  (ru) 
o  C  i  r 


(10) 


The  x,y-momentLm  equations  for  both  the  gas  and  liquid  phases 
are  given  by 
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where  S'=  - ^  Dv/mk  <uq  ‘  UW5'  (  °k  the  dra9  force  term 

defined  below),  is  thedropiet  momentum  coupling  term.  G  is  the 
body  force  and  o—  is  the  cylindrical  viscous  stress  given  by 


o  -(2-^)J.Vr 


is  the  turbulent  viscosity 


Pj-i-fTjAMciefut 


where  m  is  the  mass  of  droplets  evaporated  or  coalesced  or 
absorbedw?to  the  liquid  surface. 

Momentum  conservation  is  given  by: 


( J- Uk)  <■  Vr * Gk 
dc  mk  r 


where|def  u|is  the  magnitude  of  the  deformation  tensor . 

For  axisymmetric  flow  with  swirl,  the  two'dimensional 
Navier-Stokes  equations  are  supplemented  with  the  angular 
momentum  equation 

6  (rptv)  „  1  y  (ir^pwii)  -  —  V.  (r?)  *S% 

At  r  1 


where  is  the  body  force  or  buoyancy  term.  Uhen  a  swirl  component 
is  present,  it  is  given  by 


dW„  D„  Wk 

- - -  (w-Wk)  -Uk  — 

dc  m.  '  *  r. 


?-prJV(-p)  (15) 


Finally  the  equation  for  the  internal  energy  is  given  by: 


while  the  droplet  energy  conservation  is  given  by 


dHk  ,  ,  dink 


HsILtlv.  ( rp’Xu)  -  -  —  V.  (ru)  *? .  V(  iv/r )  ♦  —  u.  Vr 


1 / r  7 ■ l r“)  *  Q 


where  the  heat  flux  consists  of  both  conduction  and  diffusion 


where  L^  is  the  latent  heat  of  vaporization  and  qk  is  the  rate  of 
heat  transfer  to  the  droplet  from  the  surrounding  gas. 

19 

Following  Duckowicz  ,  a  Monte  Carlo  droplet  group 
representation  is  used,  i  .e. ,  each  droplet  parcel  actual  ly  represents 
a  certain  number  (can  be  fractional)  of  particles  of  identical 
characteristics  in  radius,  velocity,  temperature  etc.  Except  in  the 
case  of  very  few  parcels,  this  spray  is  expected  to  statistically 
behave  like  the  actual  droplet  cloud. 

These  equations  for  the  droptet  phase  combined  with  those  given 
earlier  for  the  continuum  gas,  l iquid  phases  describe  the  flow  ful ly 
once  the  interaction  terms  of  mass,  momentum  and  energy  exchange 
between  the  droplets  and  gas/liquid  are  defined. 


j  -  fC?T-pgDY.hkV  (  — ) 

*  P  j 


Alt  the  source  terms  in  equations  1,10,13  and  15  refer  to  the 
production  rate  per  unit  total  volume  Vg  +  and 

<19> 


The  Liquid  Oroplet  Phase  Eouations 

The  droplets  are  tracked  in  a  lagrangian  fashion.  Thus 
their  mass  conservation  is  given  by 


LIQUID  DROPLET  FORMATION 


The  first  problem  that  will  be  tackled  is  the  modelling 
of  the  liquid  droplet  formation  at  the  gas/liquid  interface.  The 
rate  of  generation  of  the  nurtier  densities  of  various  sizes  of  liquid 
droplets  at  the  gas/liquid  boundary  is  modelled  by  combining  the 
local  mass  strippping  rate  formulation  of  Mayer 
''  ,  with  the  Sauter  mean  diameter  fonrxjla  of  Ingebo  along  with 
a  probability  distribution  function  (pdf)  size  distribution 

around  this  mean  diameter,  as  per  Dukowicz  ’  . 

Liang  etal  have  found  that  Mayer's  atomization  theory  was 

quite  successful  in  modelling  the  liquid  fuel  spray  formation  in 
rocket  fuel  injectors. 

Based  on  Mayer's  theory^7  of  capi  1 1 ary  wave  ampl  i  f  icat  ion,  this 


3 


mass  stripping  rate  Mfl  is  given  by 


<2‘> 


where  the  constant  C  was  found  to  be  between  0.2  and  0.3.  vr  is 
the  relative  velocity  between  the  gas  and  liquid  at  the 
gas/liquid  inerface. 

This  quantity  M  will  be  equal  to  the  sum  of  the  masses  of 
all  the  droplets  N  of  diameter  d  formed  at  the  boundary,  namely. 


—  - (0.198,0.759,1.0,1.23,2.30) 
d. 


(29) 


Liquid  Droplet  Drag 


For  liquid  droplets  of  molecular  size  one  may  assune  no  veloci  ty 
lag  between  the  fluid  and  droplet  phases.  It  is  in  this 
appoximation  that  one  has  the  LHF  (Locally  Homogeneous  Flow)  model 
discussed  earl ier.  However,  for  l  iquid  droplets  of  significant  size 
in  the  flow  there  will  be  a  velocity  lag  between  the  carrier  gas  and 
the  liquid  droplets.  The  effect  of  this  velocity  lag  between  the  gas 
and  liquid  droplet  is  characteri zed  by  the  momentun  coupling  term 
S’in  equation  (11)  and  the  drag  force  term  in  equations  20,  21. 
It  is  given  by 


M.-EpjJV*1 


(25) 


0J,-6it|irJt*  —  r/C0|u,| 


(30) 


here  the  drag  coefficient  CD  's 

dependance  on  the  Reynolds  number  is  given  by 


Co  1 

24  (Red)‘°-S4 

Red  <  80 

where^  is  the  density  of  the  liquid. 

0.271  (Red)0'217 

80  <  Red  <  IQ4 

(3#) 

Ingebo  (20)  has  studied  the  gas  density  effects  on  liquid 

droplet  sizes  and  found  that  the  Sauter  mean  diameter  is  given  by 

= 

2. 

Red  <  104 

the  formula 

This 

variation  in  the  drag 

coefficient  with 

the  Reynolds 

-^-4 .33x10-11  (WeRep^)  0 0,5  (26> 


approximately  accounts  for  effects  of  the  flattening  of  the  droplet 
and  the  resultant  increase  in  drag  coefficient  as  the  droplet 
deforms. 


where  Ue,  Re  are  the  Weber  and  Reynolds  numbers,  do  is  the 
injector  orifice  .diameter  and  d  is  the  Sauter  mean  diameter 
defined  by  i  ndj  ■  vm  >s  the  RMS  velocity  of  the  gas 
molecules. 

As  for  the  probability  distribution  function  for  the 
particle  sizes,  Ouckowicz^0  found 


fr(r)  «-^-exp(-^~] 

d.  d. 


(27) 


where  r=d/2  is  the  radius  of  the  liquid  droplets. 

Equations  23-26  can  be  combined  together  to  determine  the 
nurtber  density  of  a  particular  droplet  size  by  noting  that  the 
nuiber  of  droplets  dN  of  radius  r  in  the  interval  dr,  located  at 
position  x,  in  the  volune  interval  dx,  and  with  velocity  u  in 
the  interval  du  at  time  t  is  given  by 


dN~f(r,x,  up,  c)  didxdup 


(28) 


Dispersion  of  liquid  droplets  due  to  turbulence 
19 

As  discussed  by  Dukowicz  ,  the  effects  of  turbulence  or.  the 
diffusion  of  the  liquid  droplets  can  be  described  in  two  ways.  One 
way  is  to  assume  a  droplet  diffusion  equation  with  a  turbulent 
particle  diffusivity.  The  evaluation  of  this  turbulent  particle 
diffusivi  ty  requires  turbulence  model  l  ing  to  evaluate  the  properties 
of  the  turbulence  field  as  done  by  Hostafa  and  Elghobashi^  . 

As  discussed  by  Dukowicz,  an  alternate  possibility  (which  is 
adopted  in  this  paper),  is  to  attempt  to  model  the  diffusion 
directly,  as  a  random  walk  of  particles  acted  on  by  the  turbulent  gas 
velocity  field.  Assuming  isotropy  and  a  Gaussian  velocity 
distribution  for  the  gas  turbulence,  random  instantaneous  turbulent 
velocities  for  the  liquid  droplets  may  be  selected  from 


•  Y)  erf~l{ |Jf|,  |5i) 


(32) 


where  X  and  Y  are  random  variables  uniformly  selected  from  the  range 
(-1  ,  +  1).  £  is  the  turbulent  kinetic  energy  level  (usually  specified 
as  a  fraction,  say  10X,  of  the  mean  flow  kinetic  energy).  This 
velocity  perturbation  acts  for  a  duration  equal  to  T  (another 
empirical  parameter),  which  can  be  thought  of  as  the  time  it  takes 
for  a  typical  droplet  to  traverse  the  width  of  a  typical  turbulent 
eddy. 

The  turbulent  dispersion  of  the  liquid  droplets  is  thus 
dependant  on  the  empiricism  involved  in  defining which  need  to 
be  evaluated  by  correlating  the  predictions  of  this  model  to 
experimental  observations. 


However  for  computational  purposes,  this  function  is 
represented  by  five  distinct  groups,  with  one-fifth  of  the  liquid 
droplet  mass  being  assigned  in  random  order  to  each  of  them. 
These  five  ratios  of  actual  size  to  mean  drop  size  are  taken  from 
Ingebo: 


NUMERICAL  SOLUTION 

The  numerical  solution  scheme  used  for  the  solution  of  the  fluid 
and  droplets  equations  given  above,  isa  modification  and  extension 
of  the  multi-step  ICE-ALE  solver  code2  '  '  One  can  refer  to  them 
for  complete  details.  Only  a  few  salient  features  are  discussed 


4 


below. 

The  scheme  is  nominally  first  order  accurate  in  both  space 
and  time.  In  the  time  dimension,  the  scheme  is  essentially  a 
split  time-step  marching  procedure,  evaluating  alt  quantities 
explicitly  using  the  most  updated  values  of  the  f*eld  variables 
upto  that  point  in  the  solution  cycle.  The  only  exception  is 
that  of  solving  for  the  ftuid  accelerations  due  to  pressure.  In 
view  of  the  strong  coupling  between  the  pressure  and  velocities, 
the  velocity  updates  are  made  implicit  with  pressure,  thus 
overcoming  Courant  stability  restriction. 

One  of  the  important  questions  that  arises  in  analyzing 
particle- laden  flows  is  that  of  the  evaluation  of  pressure.  For 
particles  of  molecular  size  or  in  the  LHF  approximation,  the 
effect  of  these  particles  is  treated  as  though  it  is  another 
species  of  a  certain  mass  fraction,  thus  entering  into  the 
calculation  of  the  value  of  density  in  the  gas  equation  p=o  RT. 
For  particles  of  more  significant  size  where  they  are  treated  as 
a  separate  phase,  their  contribution  to  pressure  may  be  modelled 
after  Statistical  Mechanics  of  continua  if  they  are  assumed  to 
have  a  completely  random  position  and  velocity  as  the  molecules 
in  the  Classical  Statistical  Mechanics.  Mo  such  analysis  seem  to 
have  been  done  todate.  The  other  way  is  to  calculate  the  fluid 
pressure  by  excluding  the  particles  in  a  given  volume,  thereby 
reducing  the  effective  areaa  over  which  the  pressure  is  applied. 
This  is  the  method  adopted  in  our  work. 

The  major  steps  in  the  solution  procedure  are: 

Step  1.  Explicit  lagrangian  evaluation  of:  Liquid  Droplet 
formation  and  transportation;  Alt  gas-droplet  coupling  terms;  All 
viscous  diffusion  terms;  Pressure  gradient  and  body  force  updates 
of  velocity 

Step  2.  Implicit  updating  of  velocities:  Point  pressure 
iterative  correction  for  fluid  velocities;  Implicit  coupling  of 
droplet  velocities 

Step  3.  Evaluation  of  convective  terms:  Gas  &  Liquid 
mass,  energy  convection;  Momentum  convection 


RESULTS 


In  this  section,  the  results  of  some  simulations  are 
presented.  This  code  is  time-accurate  and  is  suited  for 
modelling  transient  phenomena.  The  physical  duration  of  the 
simulation  run  is  divided  into  a  large  number  of  time-steps  or 
cycles  whose  size  is  continually  being  modified  on  numerical 
stabiliity  considerations.  For  fixed  boundary  conditions,  a 
steady-state  solution  is  obtained  by  running  the  problem  long 
enough  until  all  transient  effects  of  the  flow  field  died  out. 

The  principal  interest  in  this  phase  of  the  work  is  to  see 
how  well  the  models  predict  the  fluid  properties  across  the 
mixing  zone  and  along  the  combustor  centerline,  the  flow 
boundaries  (Fig. 5),  the  particle  boundaries  and  their  velocities. 
The  results  presented  here  are  for  an  oxygen  gas  jet  injected 
into  a  liquid  hydrogen  bath.  The  injector  inner  dia  is  0.089" 
and  the  combustor  diameter  is  0.62“  and  is  12"  long.  The  gas  jet 
velocities  are  varied  from  0.6e4  cm/sec  to  9.6e4  cm/sec.  The 
results  presented  are  for  200  cycles.  The  computational  grid  i*n 
the  combustor  is  shown  in  Fig. 6. 

Figs.  7a,  7b,  7c  show  the  axial  velocity  distributions 
across  the  mixing  zone  for  three  jet  velocities,  3.57e4,  5.93e4 
and  8.11e4  cm/sec  respectively.  The  velocity  is  normalized  by 
the  maximum  axial  velocity.  Each  of  the  figures  contains  the 
distributions  for  several  axial  locations  in  the  combustor.  It 
is  seen  from  the  figures  that  the  maximum  axial  velocity 
(V/Vmax*1)  ii  reached  downstream  of  the  jet  exit  plane  where  the 
extent  of  the  constant  velocity  region  i.e.  the  unmixed  core  of 
the  jet  has  ben  reduced  from  the  maximum  (x=11)  to  (x=2)  at  an 
axial  locatio  of  about  7.  The  velocity  distributions  are 
different  for  different  jet  exit  velocities. 

Fig.  8  shows  a  typical  velocity  distribution  along  the 
combustor  centerline.  The  maximum  axial  velocity  for  this  case 
is  5.93e4  cm/sec.  The  figure  shows  that  the  jet  accelerates  upto 
about  10  axial  locations  before  peaking  up.  At  about  y=25,  the 


jet  completely  dissipates  its  momentum  and  the  jet  velor'.ty 
approaches  zero. 


Figs.  9a,  9b,  9c  show  the  distribution  of  the  liquid  fraction 
F  across  the  mixing  zone.  F  =  1 . 0  means  fully  gaseous,  while  M0.G 
means  fully  liquid  and  for  other  values  of  F  in  between,  part  of  the 
cell  at  that  location  is  liquid  and  part  of  it  is  gas.  For  the 
lowest  speed  of  the  jet  of  3.57e4  cm/sec  shown  in  Fig.  9a,  one  sees 
changes  in  F  only  towards  axial  location  of  y=9.  This  means  that  the 
gas  jet  is  not  mixing  very  much  with  the  surrounding  liquid.  As  the 
speed  of  the  gas  jet  increased,  Figs.  9b  (Vmax=5.93e4  cm/sec),  and 
Fig. 9c  (Vmax  =  8.11e4  cm/sec),  the  changes  in  F  are  seen  earlier. 

The  variation  of  the  liquid  fraction  F,  along  the  combustor 
centerline  is  shown  in  Figs. 10a,  10b.  In  Fig.  10a  with  the  lower  jet 
speed  of  2.35e4  cm/sec,  the  jet  center  remains  fully  gaseous  for  a 
longer  axial  distance  than  the  higher  speed,  5.93e4  cm/sec,  case  of 
Fig.  10b.  This  means  that  the  higher  jet  speed  is  conducive  to  more 
mixing  between  the  gas  jet  and  the  surrounding  fluid. 

Finally  the  location  of  the  l iquid  particles  within  the  mixing 
zone  are  shown  in  Fig.  1  la  (Vmax=4 ,8e4  cm/sec)  and  F ig.  1 1b  (Vmax=3.6e4 
cm/sec),  in  these  figures  x=0  corresponds  to  the  jet  or  combustor 
centerline  while  y=0.3  corresponds  to  the  combustor  face  plane.  It 
is  seen  from  the  two  figures  that  all  the  particles  are  located 
farther  away  from  the  jet  centerline  and  farther  downstream  along  the 
axis.  The  particles  closest  to  the  centerline,  x=0,  are  farthest 
along  the  jet  centerline.  In  other  words  the  particle  boundary  peaks 
away  from  the  oxygen  rich  center  region  of  the  gas  jet.  As 
discussed  earlier,  the  model  for  the  turbulent  dispersion  of  the 
liquid  droplets  contains  certain  empirical  constants,  namely,  the 
resident  time  t  and  the  turbulent  kinetic  energy  fraction.  A 
parametric  variation  of  these  variables  need  to  be  carried  out  to  get 
complete  picture.  It  is  proposed  to  do  so. 


CONCLUSIONS 

Models  are  developped  to  describe  the  formation  and 
characteristics  (number  density/diameter  distributes)  of  liquid 
droplets  at  the  gas/liquid  interface  of  a  gas  jet  injected  into  a 
liquid  bath  at  rest.  Their  dispersion  by  the  carrier  phase  (gas) 
turbulence  as  well  as  the  effects  of  velocity  lag  between  the  gas  and 
droplets  is  also  modelled. 

Predictions  of  axial  and  radial  velocity  distributions  and 
.  ,q':’d  fractions  are  found  to  be  realistic.  The  location  of  the 
l iquid  droplets  wi thin  the  mixing  zone  is  shown  to  peak  away  from  tne 
oxygen  rich  center  jet.  As  the  turbulent  dispersion  model  for  the 
liquid  droplets  contains  some  empirical  constants,  it  is  considered 
important  to  do  some  more  simulations  with  a  parametric  variation  of 
these  costants. 

Work  is  in  progress  to  correlate  the  predictions  of  this  model 
with  experimental  results  of  injection  of  gas  jets  into  liquid. 
Unfortunately  there  is  a  paucity  of  data  on  the  liquid  particle 
distributions  for  such  flows.  It  is  hoped  that  this  experimental 
work  will  be  forthcoming  in  the  near  future. 
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